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Abstract 

A special class of Dirac-Pauli equations with time-like vector potentials of 
external field is investigated. A new exactly solvable relativistic model de- 
scribing anomalous interaction of a neutral Dirac fermion with a cylindrically 
symmetric external e.m. field is presented. The related external field is a su- 
perposition of the electric field generated by a charged infinite filament and the 
magnetic field generated by a straight line current. In non-relativistic approx- 
imation the considered model is reduced to the integrable Pron'ko-Stroganov 
model. 
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1 Introduction 



Exact solutions of relativistic wave equations are both very rare and important. 
Firstly they provide explicit solutions to concrete physical problems free of inaccura- 
cies and inconveniences of the perturbation theory. Secondly, such exact solutions can 
serve as convenient basis sets for expanding the solutions of other physical problems 
which are not necessarily exactly solvable. 

A good survey of exact solutions of relativistic wave equations can be found in 
Notwithstanding this book has been published as long ago as in 1990, it continues to 
be a good information resource on exactly solvable relativistic systems for particles 
with spins and 1/2. Surely this collection is not exhausting: many new results 
have been obtained during two latest decades, including the problems for the Dirac 
equation in lower dimension space and the problems for neutral Dirac particles. 

Exact solutions of the Dirac equation describing electrically neutral particles with 
non-minimal interaction with an external electromagnetic field are noteworthy. Phys- 
ically, such solutions have a large applicative value since they can be used to model 
the motion of a neutron in realistic situations. In particular they have relations to the 
nuclear reactors security problems. Moreover, magnetic trapping of neutrons is a sub- 
ject of direct experimental studies, refer, e.g., to [2]. Mathematically, the anomalous 
interaction terms depending on tensor fields present in such equations, dramatically 
reduce the number of problems which can be solved using complete separation of 
variables. In addition, just neutral particles anomalously coupled to external mag- 
netic field give rise of the Aharonov-Casher effect [3] with its interesting physical and 
mathematical aspects. 

The very possibility to exactly solve a quantum mechanical problem stems from 
the existence of a dynamical symmetry which is more extended than the geomet- 
ric symmetry of the problem. The famous examples of such systems are the Kepler 
problems and isotropic oscillator whose dynamical symmetries are defined by groups 
SO{A) and U{3) respectively. One more well known example is the Pron'ko-Stroganov 
(PS) problem [1] which describes anomalous interaction of a non-relativistic electri- 
cally neutral particle of spin 1/2 with the field of straight line current. The related 
dynamical symmetry of negative energy states is described by group 5*0(3) while the 
geometrical symmetry of the system is reduced to the rotation group in two dimen- 
sions, i.e., to 5*0(2). Let us stress that the formulation of the PS problem uses the 
Schrodinger-Pauli equation for neutral particles, i.e., it is essentially nonrelativistic. 

Paper was followed by a number of publications devoted to exactly solvable 
problems for neutral particles. In particular, the supersymmetric aspects of the PS 
model were investigated in [5]-[7], more realistic models based on magnetic field pro- 
duced by current of thin filament were discussed in [8] and [6] following an approach 
non-relativistic. A rather completed study of the Dirac-Pauli equation for neutral 
particles can be found in [H], the case of purely electric time-independent external 
field was studied in [10]. However, an exactly solvable relativistic analogue of the PS 
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problem was not known till now. 

We can add that searching for exact solutions of Dirac equation belongs to ever- 
green problems, apparently the most recent result in this field can be found in 
Exactly solvable two-particle Dirac equations are discussed, e.g., in [12]. For exact 
solutions of relativistic wave equations for particles with higher and arbitrary spins 
see [T3j| and ^14j . 

In the present paper we discuss a certain class of relativistic problems describing 
anomalous interaction of the Dirac fermion with an external electromagnetic field. 
The considered equations admit an effective reduction to equations invariant with 
respect to the 1-1-2 dimensional Gahlei group. We show that this class includes new 
exactly solvable relativistic systems and we study in detail one of them, namely the 
one closely related to both the Pron'ko-Stroganov problem and relativistic problems 
for the neutron interacting with an external field. The corresponding external field 
is a superposition of the magnetic field generated by a straight line constant current 
and the electric field of a charged infinite filament. We will show that this system 
can be treated as an exactly solvable relativistic analogue of the Pron'ko-Stroganov 
problem. 

In spite of that our main goal is to present an exactly solvable problem for neutral 
fermions, for the sake of generality we consider also more general problems with both 
the minimal and anomalous interactions. We also indicate an exactly solvable model 
of this general type, see Section 7. 



2 Dirac-Pauli equations and reduction 
SO(l,3)^HG(l,2) 

Consider the Dirac-Pauli equation for a charged particle which interacts anomalously 
with an external electromagnetic field: 

(7%-m-A^'^'^F^,)V^ = 0. (1) 

Here tc^ = — eA^, = i-£^,A^ being components of the vector-potential of the 
external electromagnetic field, j'^ are Dirac matrices satisfying the Clifford algebra 



g'^" is the metric tensor whose nonzero elements are g'^^ = —g^^ = —g'^'^ = —g^^ = 1, 
S^" = i {j^j'' — 7^7^) /4 is the spin tensor, F^^ = n^] is the tensor of the electro- 
magnetic field such that F^a = —Fa, Fat = SahcBc where Fa and Be are components 
of vectors of the electric and magnetic field strengthes. In addition, e and A denote 
particle charge and constant of anomalous coupling. The latest is usually represented 
as 

A = giiQ (3) 
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where /iq is the Bohr magneton and g is the Lande factor. 
We use Heaviside units with h = c = 1. 

Equation ([1]) describes both the minimal and anomalous interaction of a Dirac 
fermion with an external electromagnetic field. Setting in ([1]) A = and supposing 
e 7^ we come to an equation describing the anomalous interaction only, while for 
e = 0, A 7^ we have a Dirac-Pauli equation describing a neutral fermion. 

Equation ([1]) is transparently invariant with respect to the Lorentz group SO (1,3) 
which transforms time and space variables xo,xi,X2,X3 between themselves. Among 
the subgroups of this group there is the homogeneous Galilei group HG(1,2) which 
includes the transformations of variables r, Xi,X2 where r = {xq + x^)/2 (for all non- 
equivalent subgroups of S0(l,3) and of the Poincare group see 

To search for exactly solvable problems based on equations ([T]) we restrict ourselves 
to special class of external fields, which makes it possible to expand solutions of ([T]) 
via solutions of reduced equations invariant w.r.t. group GH(1,2). In other words, 
we will discuss such external fields for which these reduced equations be integrable. 

To this end we first suppose that the vector-potential A = {Aq, Ai, A2, A^) be 
light-like, i.e., 

A^A'^ = 0. (4) 

This condition can be always satisfied up to gauge transformations and so it does not 
lead to any loss of generality. Then we restrict ourselves to the vector-potentials of 
the following special form compatible with (jl]): 

A=(<^,0,0,<^) (5) 

where if is a function of time and spatial variables. In addition we suppose that ip 
depends on three variables only, namely, 

if = ip{T,Xi,X2). (6) 

Vector-potentials ([5]), ([6]) satisfy the Lorentz gauge condition p^A^^ = identically 
and the invariants of the related external field are both equal to zero, i.e., 

F^.F'^^ = and i^^.^^F^'^F'"^ = 0. (7) 

For convenience we fix a nonstandard realization of the Dirac matrices and set 
o_ f I \ 3 f I \ . f -a^ 



^ = > / 0^'^ = ^-/ ^ =\ 

where a = 1,2, are Pauli matrices, and / are the 2x2 zero and unit matrix 
correspondingly. Then 

^ - 2 I ' ^ - 2 U. • 
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System ([T]) with a particular class of the vector-potential given by equations 
and ([HD is homogeneous with respect to the difference of independent variables Xq — x^. 
Thus it is convenient to rewrite it in the light cone variables 

'T = 2^X0 + X3) and ^ = xq- X3. 
As a result we obtain: 



= (7^7f^ - m - Xr]^F^) ^ = (10) 
where = a = 1,2, 

70 = 70 + 73, 73 = lilO- 73), 7a = 7a, ^a = K7o7a + 737a), 

(11) 

vTo = i|: - 2ev9, TTs = P3 = -i^, vr^ = = -i^ 

and summation w.r.t. repeated indices fi and a is imposed over the values /i = 0, 1, 2, 3 
and a = 1,2 respectively. In addition, we impose on solutions of (fTOj) the standard 
condition of square integrability and ask for ip ^ when Xa 0. 

Operator P3 commutes with L and so is a constant of the motion for equation 
( TTP]) . Let us expand solutions of this equation via eigenvectors ipM of P3: 

Psi^M = MijM Vm = exp(iM0^(r,xi,X2). (12) 
Let us denote 

V'(r,xi,a;2) = ) ( (13) 

where p and x are two-component spinors. Substituting f|T2l) and f|T3|) into f fTOj) and 
using realization of 7-matrices we obtain the following system: 

{ia^Pa -m)p+ (i-^ - 2eip - Acr„F„ J x = 0, (14) 



2Mp-{i(T^p^ + m)x = 0. (15) 

It is easy to convince oneself from eqs. fll4p . (1151) that, without loss of generality 
we can assume that M cannot take the zero value. Indeed, setting M = in f|T5l) we 
reduce it to the equation for x which does not have non-trivial normalizable solutions. 
Then, equating x to zero in f|T^ . we obtain the equation for p whose normalizable 
solutions are trivial also. 

It is interesting to note that this system is nothing but a (l+2)-dimensional version 
of the Galilei-invariant Levi-Leblond equation [16j with anomalous interaction, as can 
be immediately deduced by comparing with equation (52) for e = A; = in [T7] . 
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Galilei invariance of system iHM . ffT5]) can be easily verified since the Galilei boost 
for coordinates Xa- 

where Va are the boost parameters, keeps this system invariant provided Fa and ipM 
cotransform as 



where C is an arbitrary constant. Transformations f|T6l) belong to a representation of 
the Galilei group in (l+2)-dimensional space. 

Thus the Dirac-Pauli equation ([T]) can be reduced to the set of Galilei invariant 
equations ffUl) . ffT^ provided the related vector-potential of the external electromag- 
netic field belongs to the class described by formulae ([5]) and ([6]). Solving eq. (|T5i) 
for p under the condition M 7^ and substituting it into eq. (|T^ we obtain the 
Schrodinger-Pauli equation for the two-component spinor x'- 



where ^0 = ^ and = p\ + pi- 

Surely the non-relativistic equations (ITTjl are more easy to handle then the initial 
relativistic equation ([1]), since they include smaller numbers of dependent and inde- 
pendent variables. In particular, a number of exactly (and quasi-exactly) solvable 
Scrodinger-Pauli equations (ITTI) is well studied, and many of them can be used to 
construct solvable relativistic problems using the scheme inverse to the previously 
proposed. 

In Section 4 we use this idea to generate a relativistic analogue of the PS problem. 

3 Cylindrically symmetric potentials 

Consider in more details a physically interesting subclass of equations ([I]), ([5]), ([6]) 
when the corresponding potential (p depends on the square x'^ = x\ + x\ of 2-vector 
X = (xi,X2) and is independent on r. The related reduced equation f|T71) takes the 
form 



F, 



a 




(16) 




(17) 




(18) 



where 



H = eo + 




+ 2eip + X 



(TaXa 

X dx 



(19) 
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Equation f[T5]) has three additional constants of motion, namely, 



Po = i 




(20) 



where Ri is the reflection operator which acts on x as follows: 
RiXir, xi, X2) = xi'T, -xi, X2). 

Operators Pq and J12 are generator of shifts w.r.t. variable r and rotation gener- 
ator correspondingly. They commute with the Hamiltonian f|T9l) and between them- 
selves. Expanding solutions of f|T8|) via complete sets of eigenfunctions of Pq and Ju 
it is possible to separate variables in this equation. 

Operator Q represents a discrete symmetry w.r.t. the reflection of the first coor- 
dinate axis. It commutes with Pq and H (|T9|) but anticommutes with J12. It follows 
from the above that eigenvalues of Pq and H should be degenerated w.r.t. the sign 
of eigenvalues of Ju- 

Notice that equation (1181) admits other discrete symmetries like reflections of X2 or 
both Xi and X2- But all such additional symmetries are either rotation transformations 
or products of reflection Q and rotations. 

Let us separate variables in equation f|T8|) . First we define the eigenvectors of Pq 
which have the following form: 



where H is Hamiltonian f|T9|) . 

In addition to the coupling constants e and A, equation fl221) includes two param- 
eters, i.e., 6 and M. For a fixed non-zero M this equation defines an eigenvalue 
problem for e. 

Now we can use the symmetry of (122|) w.r.t. the rotation group (whose generator 
is J12) to separate radial and angular variables. To do this we rewrite equation (122|) in 
terms of angular variables, i.e., set Xi = xcosO, X2 = xsin6, r = 2MXx (where A is a 
normalizing parameter), and expand x via eigenfunctions of the angular momentum 
operator J12: 



where Ck are constants, (pi and 02 are functions of r and summation is imposed over 
the repeated indices k = 0, ±1, ±2, • ■ ■ . 



Xe = exp[-ier)x[^). 



(21) 



Then, substituting eq. fl2Tl) into eq. ffTSl) we obtain the equation 



^X = Hx 



(22) 




(23) 
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In the following we restrict ourselves to solutions Xk which correspond to non- 
negative values of k. Then solutions with k negative would be obtained by acting on 
Xk by operator Q fl20l) . 

Substituting eq. fl23l) into eq. fl22l) we come to the following system: 

H,<f> ^ + k{k - aa)^ + 2e^ + a.j^'^ <f> = e<j> (24) 

with (f) = column(</)i, (1)2), and 

e={e- eo)/2M~X\ (25) 

Thus we reduce f l22l) to the system of two ordinary differential equations for radial 
functions, given by formula fl2^ . Its solutions must be normalizable and vanish 
at r = 0. For some types of potential (f (and particular restrictions imposed on the 
coupling constants e and A) this system is integrable and its solutions can be expressed 
via special functions. In the following section we consider an example of integrable 
equation fl24p which corresponds to a neutral particle interacting anomalously with 
an external field. 



4 Relativistic analog of PS problem 

Let us set e = in (1141) . (1151) and choose the following particular realization for the 
potential (p: 

(f = u log(x) (26) 

where is a constant. Then the related equations (ITSI) and are reduced to the 
following forms 

^X= {wT7 + ^ — e =e-eo (27) 



2M x2 



and 



Hk^= i-— + k{k-a3)^ + (Ti-]^ = i(p (28) 

correspondingly, provided we set A = cuA. 

The electromagnetic field whose potential is defined by relations and (|2^ has 
a transparent physical meaning. Namely, it is a superposition of the electric field 
E = {El, E2, E3) whose components are 

Ei = u^, E2=u;^, Es = (29) 
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and the magnetic field B = {Bi, B2, B2) with 



5i = -cu^, ^2 = 0^^, 53 = 0. (30) 

Such an electric field can be identified as the field of a charged infinite filament 
coinciding with the third coordinate axis. Let us designate the charge density of this 
filament by p then the coupling constant uj should be equal to 2p. On the other hand, 
the magnetic field B is nothing but the field of a straight line constant current J 
directed along the third coordinate axis provided the coupling constant uj be equal 
to 2 J. Of course the related charge density and current should be equal between 
themselves, i.e., 

j=p = uj/2. (31) 



Let us show that equation fl27j) is exactly solvable and find its solutions. The 
simplest way to prove integrability of fl27|l is to make the unitary transformation 

,_,,.,„,_(|.,,-.).,(,_(|l,X^))m (32, 

where U = "^(1 ~ i^s)- As a result we reduce (127|) to the following form: 

( r,\ SiX2 - S2Xi \ 

(,2M-2A Jx (33) 



e'x 



where 5*1 = |cri and 5*2 = ^(12 are fermionic spin matrices and in accordance with 
eqs.©, ([26D and ([31]) A = Acj = 2gpoj- 

For a fixed M and up to the value of the coupling constant equation fl33l) coincides 
with the Schrodinger equation for a neutral particle minimally interacting with the 
field generated by an infinite thin current filament (in our case the standard coupling 
constant is multiplied by factor 2). This equation was studied in numerous papers 
starting with [4j and continuing with [5]-[7] and many others. It has the following 
nice properties: 

• equation (133|) admits a hidden dynamical symmetry with respect to group 5*0(3) 
for negative eigenvalues e, group S0{1, 2) for e positive and group E{2) for e = 

• it possesses a hidden supersymmetry [5]; 

• using any of the above mentioned properties the equation can be integrated in 
closed form [l]-[6]. 

Since equation (1271) is unitarly equivalent to eq. (l33i) it has the same properties. 
In particular, eigenvalues e' are the same in both equations fl27j) and (l33l) . 
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5 Relativistic and quasi relativistic energy levels 



In the next section we will present exact solutions of equation fl27|) for coupled states 
and define the related eigenvalues e' . In fact these eigenvalues are well known, and 
using directly results of paper [4] (or of the papers [5]- [8]) we can immediately write 
e' in the following form: 

■ 2^'^*^ (34) 



Ar2 

where iV is a positive natural number. 

Eigenvalues (IMIl are degenerated since they do not depend on eigenvalues k of the 
angular momentum operator Jyi- The degeneration factor is equal to 2^ + 1, and the 
quantum number can be represented as 

AT = 2(n + A;) + 1 (35) 

where n is a natural number (1], [6]. 

Using flM|) we already can find energy levels for the initial relativistic problem. 
Indeed, since -Pq = Po + and P3 = |(po ~ Ps); it is possible to write analogous 
relations for eigenvalues of po? of Ps aiid e, M: 

e = E + K, 2M = E- K. (36) 

Then, using definitions (1271) and for e' , E and M we find from the relativistic 
energy spectrum: 



E=-^\^ + \l + '^' + ^\ (37) 




where R = K/m and A = 2fiogj. 

We see that in spite of the fact that the neutron motion along the third coordinate 
axis is free, the third component of momentum k makes a rather non-trivial contri- 
bution into the values of energy levels (jSTj). In accordance with flTTI) E > n, and so 
the condition M 7^ in fact is satisfied. 

The most simple expression for energy levels corresponds to the particular value 
K = 0: 

E = (38) 
1 + * 

which for small A becomes 

' 2A^ +'"""' A^^ + ■ ■ ■ ■ ^ ^ 
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Up to the rest energy term m the approximate energy levels fl5I?l) are exactly 
the same as in the non-relativistic PS problem [3], In particular both the 

approximate and exact levels given by equations (137|) . (l38l) and (l39l) are degenerated 
with respect to eigenvalues k of the third component of angular momentum which is 
a constant of motion for the considered system. Like in [3] this degeneration is caused 
by a hidden dynamical symmetry of the system. 

Both A and k can be treated as small parameters. Expanding E fl39p in power 
series of A and k we obtain the quasirelativistic approximation for the energy levels: 

E ^ m + 1- 40 

where 

A. = (1 - ^)A. (41) 

The first tree terms in fl40p represent respectively the rest energy, the kinetic en- 
ergy of the free motion along the third coordinate axis and the relativistic correction 
to this energy. The last (dynamical) term in fj^Ol) is quite similar to the correspond- 
ing non-relativistic term (compare with fl39l) ). but includes the corrected coupling 
constant A^ instead of A. 

Notice that since the electromagnetic field defined by relations ^ and (jH]) has no 
components in x-^ direction, the motion of the particle described by equation (ITUl) in 
this direction is free. This free motion can be quantized by imposing the periodic 
boundary condition. Then 

ItiN ~ , , 

K = ——, N = 0, ±1, ±2, • • • (42) 
Ij 

and energy levels (|57|) - fHni) be labeled by the pairs of quantum numbers N and N. 

6 Exact solutions for bound states 

To find the solutions of eq. fl27|) we use the fact that the Hamiltonian 

/ 92 1 1\ 

H.= [-^ + k{k-a,)- + a,-j (43) 

can be factorized as 

Hk = a+ak + Ck (44) 



where 



d d 1 

flfc = — + Wk, a+ = -— + Wk, Ck 



dr ^' ^ dr " {2k + lf 
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and is a matrix superpotential 

W, = -0-3 - -^ai - (45) 
" 2r 2k + l r ^ ' 

It can be verified by direct calculation that 

(9^ 1-1 
Hk = dkat + Ck = + ik + l){k + 1 - as)— + Aai- 

i.e., the superpartner Hamiltonian if^^ for Hk is equal to Hk+i- Thus the eigenvalue 
problem fl24|) possesses a supersymmetry with shape invariance and so it can be solved 
using the standard technique of the supersymmetric quantum mechanics [18]. We will 
not reproduce the related routine calculations whose details can be found in [B], [7] 
but restrict ourselves to the presentation of the solutions of eq. fl27|) . 

The ground state solutions 0(0, k; r) = column (0i(O, k; r), 02(0, k; r)) are square 
integrable and normalizable solutions of equation afc0(O, k; r) = 0. They can be 
expressed in the following form: 

^,(0,k;r) = r'+'K,{^^), ^,(0^ k;r) = ~r'+'Koi^,) (46) 

where Kq and Ki are modified Bessel functions. The corresponding eigenvalue ik in 
(127]) and ([2l]) is equal to -j2kTW 

Solutions corresponding to the first excited state, i.e., to n = 1 are 0(l,A;;r) = 
a^0(O, k + 1), or, being written componentwice 



0i(l, fc; r) = - (f + ^) 0i(O, k+l;r)- j^^M^, k + l;r) 
= i2k+i){2k+3) ^^^^-^0 ( 2Fr3 ) - (2A; + l)r''+^Ki{^^), 

02(1, A;; r) = - (f + ^) 02(0, k+l;r)- j^M^, k+l;r) 
= {2k + ?>)r'''^^KQ{,^) - (2fe+f)|2fc+3) r'^^^Ki ( 2^ ) . 



(47) 



The corresponding eigenvalue ik is equal to —j2(k^kj+Ty ~ — (2^+3)2 • 

Finally, solutions which correspond to an arbitrary value of the quantum number 
n > can be represented as 

(f){n,k;r) = a^a^^^ ■ • ■ ajl+„_i0(O, j + n; r), n = 1,2,--- 

which gives rise to the recurrence relations 

01 (n, k; r) = -|:0i(n - 1, A; + n - 1; r) - ^±^0i(n - 1, + n - 1; r) 
+ 2(k+nyi Mn,k + n-l-r), 

(48) 

(n, fc; r) = -|:02(n, k + n-l]r)- ^^±^^02(n, /c + n - 1; r) 



+ 2(fc+l)-i '^i(^^^ + ^-l;0- 
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The related eigenvalue ik is given by relations fl51|) and fl5^ . 

It is now possible to present exact solutions of the initial Dirac-Pauli equation 
defined by relations ([T]), ([5]) and (126|) . In accordance with the above such solutions 
are labeled by the main quantum number which can be expressed by equation 
( l35|) . and by eigenvalues k and k of the third component of momenta and total 
orbital momentum. Using equations ([121), (US]), (USD, (EH), (ESD and (jlSD-dlHD we find 
these solutions in the following form: 



ipn,K,k = exp(-i(Exo - 



V2nLr 



( ex^{i{k - l)9)r]i{n,k;r) \ 
exp{i{k + \)9)7]2{n, k; r) 
exp{i{k — i)^)0i(n, k; r) 
\ exp{i{k + ^)e)(j)2{n,k;r) j 



Here r = [E — x\ + E and k are given by equations (1571) and 
and 4'2{n, k; r) are functions defined by recurrence relations fH6|) . fHHI 

r/i(n, fc; r) = A (f + fj) 02(n, A;; r) + A;; r), 

r]2{n, /c; r) = A (f - ^) 0i(n, fc; r) + ;^02(n, fc; r). 



(49) 



(50) 



Solutions (jini) are normalizable and tend to zero with r ^ 0. They are defined for 
non-negative eigenvalues k of the total angular momentum while solutions for k nega- 
tive can be obtained acting on fH9|) by the reflection operator Q = ^707273-^1 where R 
is the reflection of the flrst spatial variable, i.e., Riip{xo, Xi, X2, x^) = ip{xo, —Xi,X2, x^) 
and i?i^/'(xo, r, 6^, X3) = ip{xo, ,r, —9,X3). Operator Q transforms solutions into solu- 
tions and changes the sign of k. It is reduced to operator Q (1201) when acting on the 
third and fourth components of bispinor (H9|) . 

Let us present components of solutions (H9l) for n = and n = 1 in explicit form. 
If n = then the related functions 0i(n, fc;r) = 0i(O, fc;r) and 02(0, A;;r) are given 
by equation (H^ while rii{0, k; r) and r]2{0, k; r) have the following form: 



Vi{0, k;r) - 
r]2{0,k;r) = - 



E-K 



2k+l 



+ 



'^'K,{^)-\{2k + l)r'K,{^) 
r'^'Ko (^) . 



If = 1 then the corresponding functions (f)i{l,k]r) and 02(1, fc;r) are given by 
equations fH7j) . and 



A 



rii{l,k;r) = ( A(2A; + 1)(2A; + 3)r^ + 



6fc+7 
2A:+3 



(2A; + 1) 



2fe+3 
^1 (ra) ' 



+ 



4(fc+l) 



E-K I (2A,-+3)(2fe+l) 



fc+2 



^0 (iFTs) 



^2 



;i,/c;r) 



I 2 ^ 



' 3k±A 

^2k+l ' E 
4(A;+1) 
(2fc+3)(2fc+l) 



+ ^^i2k + 3))r^^^K4^) 
r'^'K, (^) . 
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7 Discussion 



In Sections 2 and 3 we study a rather extended class of Dirac-Pauli systems which 
can be effectively reduced to a set of Schrodinger-Pauli equations. But the main 
scope of the present paper was to find an integrable relativistic formulation of the 
non-relativistic PS problem This goal cannot be achieved by a straightforward 
relativization of the PS problem, i.e., keeping the external field as being generated by 
a filament current j and changing the Schrodinger-Pauli equation by the Dirac-Pauli 
one. The simple reason is that proceeding in this way one obtains a model which is 
non integrable. 

In this paper we succeed in obtaining an integrable relativistic analogue of the PS 
model. To this end we add an extra external field generated by an infinite charged 
line. Moreover the charge density p is related to the current j by equation fl3T|) or, in 
SI units: 

P = j/c' (51) 

where c denotes the speed of light. In this way we come to the relativistic model 
defined by relations ([T]), and fl2^ . This model is indeed integrable, its exact 
solutions for bound states are given in Section 6 and the corresponding energy levels 
are discussed in Section 5. 

In accordance with (ISTj) the required charge density is very small and tends to 
zero in the non-relativistic approximation when c —* oo. Thus we claim that the 
non-relativistic limit of the model defined by equations ([T]), and (12^ is exactly 
the PS model. This property can be directly proved using the Foldy-Wouthuysen 
transformation [19j. We shall not present such a proof here but restrict ourselves 
to comparing the energy levels predicted by our relativistic and non-relativistic PS 
models. 

In accordance with fl38p in the non-relativistic approximation these levels simple 
coincide (of course up to the rest energy term), which confirms our supposition. We 
also find corrections for the PS levels in the quasi-relativistic approximation (140|) . 
Albeit the motion along the third Cartesian coordinate is free, the third momentum 
component k makes indeed a contribution into the effective coupling constant A^. 
The origin of this contribution is the anomalous interaction of moving neutral Dirac 
fermion with the magnetic field caused by the charged line. In the rest frame the 
motion along this line is effectively changed by the current which flows in the line in 
the opposite direction, which is in perfect accordance with equation fj^B . 

The contribution of k into the effective coupling constant fHTl) is small, namely, 
in SI units it is proportional to the inverse speed of light. However, it affects the energy 
levels f HU]) much more than the relativistic correction to the kinetic energy —n'^/Sm^ 
which is proportional to the squared inverse speed of light. 

Considering the exact solutions presented in Section 6 it is easy to show that 
the probability density corresponding to wave functions (l49l) admits the following 
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approximate representation: 




where the dots denote terms of order In other words, up to the terms propor- 

tional to the inverse squared speed of hght, the probabihty density looks as in the 
non-relativistic PS problem. Thus many of conclusions concerning the probability 
distributions for the non-relativistic PS problem (such as the long lifetime for bound 
states with sufficiently large A;) made in |2] and [B] can be generalized to the relativistic 
case. 

In conclusion we have presented a new exactly solvable problem for the Dirac-Pauli 
equation describing a neutral particle which interacts anomalously with a rather par- 
ticular external field given by equations fl29|) - fl3T]) having however a clear physical 
meaning. This type of anomalous interaction is the key to expand solutions of the 
problem via solutions of the (l-|-2)-dimensional Levi-Leblond equation invariant with 
respect to Galilei group. Moreover, the considered problem possesses a hidden sym- 
metry and supersymmetry which cause the (2n+l)-fold degeneration of the energy 
levels given by equations (137|) and (135|) . 

A natural question arises whether the considered relativistic problem with its 
symmetries is unique or there are other problems which can be effectively solved 
using reduction S0(l,3)— >IIG(1,2). In Sections 2 and 3 we study a certain class 
of such problems which can be effectively reduced to radial equation flMl) which is 
exactly solvable when e = and = log(a;). We believe that there are other exactly 
solvable equations ^UM and at least two of them can be immediately written down 
if we set = a/x and consider the alternative cases e = and e 7^ 0. The related 
Dirac-Pauli equations ([I]), ([5]), ([6]) can be solved explicitly. We plane to study these 
and probably other integrable models in future. 
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